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BUCKLING OF A PRESSURIZED TOROIDAL RING 
UNDER UNIFORM EXTERNAL LOADING 

By George E. Weeks 
Langley Research Center 

SUMMARY 

A theoretical investigation of the buckling of a pressurized toroidal ring under a 
uniformly distributed line load is carried out, with finite shear stiffness, extensional 
stiffness, and appropriate contributions from the internal pressure taken into account. 

The equations governing both in-plane and out-of -plane buckling behavior of the ring, 
with direction of loading included, are derived and solved. The solutions to the buckling 
equations are presented in terms of five general parameters which are functions of the 
internal pressure, external load, and ring geometry. 

The major results of the investigation are: (a) The magnitude of the buckling load 
is very sensitive to the direction of loading during the buckling deformation and to the 
ratio of cross-sectional radius to the radius of the ring; (b) the shear stiffness of the 
ring wall material and internal pressure both contribute to the effective transverse shear 
stiffness of the ring, which strongly influences the buckling loads. 

INTRODUCTION 

Expandable structures capable of being packaged into a small volume and erected 
into load-carrying structures by inflation have been used for a number of space missions. 
One such structure is an inflatable toroidal ring which has potential use as a stiffening 
element of deployable aerospace structures; for example, it may be efficient as a stiff- 
ening ring of a deployable planetary entry decelerator. 

To use this ring in a practical structure, its buckling strength under compressive 
loads must be determined. Investigations of the stability of rings under uniform loads 
may be found in references 1, 2, and 3. These studies, however, do not account for finite 
transverse shear stiffness including effects of internal pressure, which can be expected 
to be important for inflatable rings. 

The purpose of the present paper is to determine the buckling load of a pressurized 
toroidal ring subjected to a uniformly distributed radial load, taking transverse shear 
deformation into account. Since the investigations of references 1 to 3 have shown that 


the direction of load during deformation can have a significant effect on the magnitude of 
the resulting buckling load, three directions of load during the buckling process are con- 
sidered: (a) The load remains in the same direction, (b) the load remains directed toward 
the initial center of curvature, (c) the load remains normal to the deformed surface at 
every point. Both in-plane and out-of-plane ring behavior are included. 

SYMBOLS 


A = 27rrt 

A = 7ir 2 

A* cross-sectional area of toroidal membrane perpendicular to the centroidal 

axis after deformation (see fig. 3) 



e a local circumferential extensional strain in membrane 

e a g local shear strain in membrane 

E Young’s modulus 

G shear modulus 

1 area moment of inertia of membrane 

J area polar moment of inertia of membrane 

K, K determinants defined by equations (44a) and (50a), respectively 

Mj twisting moment about the centroidal axis (see fig. 1) 

n integer 

local shear-stress resultant 
p internal pressure 

P internal-pressure stiffness parameter, pAR 2 /EI 

2 



q 


external uniform line load 


q external load parameter, qR^/ El 

Q, Q determinants defined by equations (44c), (44d), (44e), (50c), and (50d) 

r radius of cross section of toroidal membrane 

R radius from origin to centroidal axis of toroidal membrane 

R radius parameter defined by equation (11) 

R radius parameter defined by equation (36) 

s* length of centroidal axis after deformation 

s* length of outer surface of deformed ring in contact with line load 

S shear -stiffness parameter, AGR^/2EI 

t membrane wall thickness 

T extensional- stiffness parameter, EAR^/eI 

uj., U 2 , ug tangential, radial, and out-of-plane rigid-body translations of cross section 
of ring 

u, v, w tangential and normal displacements of a point in the surface of the membrane 
(see fig. 1) 

Vg, Vg shear forces in the y and z directions, respectively (see fig. 1) 

W work of external load 

x, y, z rectangular Cartesian coordinates (see fig. 1) 

Z, Z determinants defined by equations (44b) and (50b), respectively 

a circumferential angular coordinate of centroidal axis 
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angle defined by equation (14) (see fig. 2) 


0 

/3 angle defined by equation (37) (see fig. 7) 

yg > y 3 shear strain of cross section in y and z directions, respectively 
r twisting- stiffness parameter, JG/EI 

AR radial distance of point on deformed outer surface from corresponding point 

on undeformed outer surface, under the line load 

AV change in enclosed volume of membrane due to deformation 

Af2 net area enclosed by deformed and undeformed membrane at its outer surface 

9 angular coordinate of ring cross section 

111 change in strain energy of deformed membrane 

11 2 change in potential energy of internal pressure 

11 3 additional energy of the shell due to buckling 

u>i, o> 2 , W 3 rotation of ring cross section in yz, xz, and xy planes, respectively 
Subscripts : 

a constant-direction loading 

b radial loading 

c hydrostatic loading 

cr critical 


ANALYSIS 

A segment of a pressurized toroidal ring configuration and the coordinate systems 
are shown in figure 1. The radius of the ring is R and the radius of the cross section 
is r. The angular coordinates a and Q are in the plane of the ring and in the plane 
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of the ring cross section, respectively. They thus define the position of a point on the 
surface of the ring. Displacements u, v, and w of points on the surface of the ring 
are related to rigid-body translations of and rotations about the centroidal axis of the ring. 

Translations u^ of the centroidal 
axis are referred to a local Carte- 
sian coordinate system with its ori- 
gin along the centroidal axis: uj 
is tangential to the centroidal axis; 

U2 and U3 lie in the plane of the 
ring cross section. Rigid-body 
rotations wj are referred to this 
same local coordinate system. 

These directions are shown in 
figure 2. 

Certain simplifying assump- 
tions have been made in this analy- 
sis: (1) The ring is considered to 
be a membrane in that the local 
bending stiffness of the shell walls 
is neglected, (2) the ring buckles 
without localized wrinkling of the 
membrane walls, (3) any cross sec- 
tion of the ring remains rigid in its 
own plane so that the cross section remains circular and the resulting deformations can 
be characterized by the six rigid-body motions of translation and rotation of the cross 
section, (4) the rigid-body rotations are small enough for the displacements on the surface 
of the toroidal ring to be represented by the vector sums of the displacements due to 
translations and rotations of a cross section. 

Derivation of Governing Equations 

The buckling equations of the pressurized toroidal ring deforming about some pre- 
stressed equilibrium configuration can be obtained from the variational equation 

611! + 6n 2 + 6II3 - 6W = 0 (1) 

where IIj is the change in strain energy, II2 is the potential- energy change of the 
internal pressure, II3 is the potential energy of the prebuckling membrane stress acting 
through the buckling strain, and W is the work done by the externally applied loads 
during buckling. 



Figure 1.- Ring-segment configuration and coordinate systems. 
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Local buckling displacements in terms of rigid-body motions of a cross section. - 
The local displacements u, v, and w of a circular toroidal ring can be expressed in 
terms of the rigid-body motions uj and aq of a cross section. These are identical to 
those given in reference 4 for a circular cylindrical beam. They are 


u = uj + rco2 cos 9 - ro >3 sin 9 

v = -rcoj + U2 cos 9 - U3 sin 9 + " w 3^) s i n 0 cos 9 

w = U 2 sin 9 + U 3 cos 9 - —h> + o^cos^# + o^sin^#) 

2 V 


(2a) 

(2b) 

(2c) 


Note that ui and are functions of a only. 

Strain ene rgy of the toroidal ring.- If it is assumed that the strain eg in the 9 
direction is negligible, the strain energy IIj for the toroidal membrane ring becomes 


n l 


t 

2 



2 


+ 


Qp 2 

Ue Q!6> 


; )(R + r sin 0)r da d0 


(3) 


The strain displacement relations for a toroidal membrane are (ref. 5), after some 
manipulation, 


u,„ + v cos 9 + w sin 9 


e a - 


_ ’O' 


R + r sin 9 


2(R + r sin 0)2 \ V ’ a 




,2 + u 2 _ 2UW,,-,, sin 9 


2uv, a cos 


(4) 


u,g v, a - u cos 9 
01 ® r + R + r sin 9 


(5) 


where commas denote differentiation with respect to the subscripts. Nonlinear terms 
are retained in e a in order to take into account the nonlinear stretching of the ring 
during buckling. However, for the theory considered here, it is adequate to include only 
the linear terms in the shear-strain displacement relation. 

These strains can be expressed in terms of the rigid-body rotations and translations 
of a cross section of the ring by substituting equations (2) into equations (4) and (5). The 
results are 
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-a 


( Clio 2 a>-* 2\ 

= 1 ui ' + U2 + r cos 0(^2’ " w i) “ r sin ®l w 3 f + "T - + o 

R + r sin e\_ x ' \ * * /_ 


2 (R 


o P 2 ( w l') 2 + ( u 3') 2 + u l 2 + u 2 , ( u 2* ” 2u l) 

+ r sin 0) z L x ' 


+ 2r sin ©(coi’ug’ + - w 3 u l) " 2r cos 0(U2' - u^coi' + o>2) 


+ r 2 o; 3 2sin 2 0 + r 2 ^ cos 2 0(w 2 + 2a>j'^ - 2r2a> 3 sin 0 cos 0(a> 2 + <*>1/) 


(6a) 


'Q!0 


= — 7 -rw/ - (R + r sin 0 )fu><> sin 0 + a)o cos 0) 

R + r sin 0 L 1 ^2 6 J 


+ sin 0(-U3*) + cos 0^u 2 ' - u-^ - r<u 2 cos20 + r sin 0 cos 0(0)3 j| 


(6b) 


where primes denote differentiation with respect to a. If equations (6) are substituted 
into equation ( 3 ) and the result is integrated over 0, the variation of the strain energy 
IIj becomes (retaining all terms up to order (r/R) 2 J; 



where 

I = 7 !r 2 t J = 27 rr 2 t A = 27 rrt 

It should be noted that the equations here are applicable to an orthotropic material since 
e^ was assumed to be zero. Therefore, to define the stiffness characteristics of the 
ring, it was only necessary to define the four effective stiffness quantities EA, El, GA, 
and GJ. 
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Potential energy of the pressurizing gas.- The work done by the internal pressure is 

the product of the pressure and the change in enclosed volume due to deformation. The 

internal pressure is assumed to remain 

v constant during buckling. Since the cross 

segment of deformed section of the ring is assumed to remain 
centroidal axis . , , , . . , , , 

circular and to move as a rigid body 

Is* 

during deformation, the change in volume 
is 

AV = \ A*ds* - 2irRA (8) 

Undeformed S* 

centroidal axis 

' x where ds * is a differential length along 
the centroidal axis of the ring after buck- 
ling (see fig. 2), A* is the area normal 
L to the deformed centroidal axis, and 

Figure 2 - Geometry of deformation of centroidal axis of ring. — „ 

A = nr*. 



Figure 2 - Geometry of deformation of centroidal axis of ring. — 


The radial, tangential, and out-of-plane displacements of the centroidal axis of the 
ring in terms of the rigid-body displacements are U2, uj, and -U3, respectively. 

The length ds * can be expressed in the fundamental form 

(ds*) 2 = dx^ + dy2 + dz2 (9 


where 


x = R sin(o! + /3) 
y = R cos(q! + /3) 


z = -u 3 


R = \j(R + u 2 ) 2 + H 2 

Substituting equations (10) and (11) into equation (9) gives 

(ds*) 2 = jl + 2j8’ + (^^R 2 ^ 2 + (di*) 2 + (du 3 ) 2 


which can be written as 


(ds *) 2 = 1 + 2/3' + (j3') 2 + ^—5- + - 2 o H 2 do<2 

R R _ 
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where 


R = r(i + ^ + ^ 5 - + . . .) 
V R 2R 2 ) 

(13) 

= tan - 1 =r—^ — = — — — + . . . 
R + U 2 R + U 2 

(14) 


Expanding equation ( 12 ) and retaining up to second-order terms in the displacements gives 


ds* = jl + Ul ^ U - + ^ 2 ^ 1 2 + ( u 2 ’) 2 + ( u 3’ f ~ 2 u 1 u 2 ^| r da ( 15 ) 


To obtain the area A* normal to the deformed centroidal axis, note that the cross- 
sectional area A is unchanged but is no longer perpendicular to the deformed axis at 

: coordinates (r, a+/3, -U 3 ) because of the shear deforma- 
tions of the ring. An illustration of this effect for only 
the shear deformation in the plane of the ring y 2 is 
shown in figure 3 . 

From figure 3 , the area A* can be expressed 
in terms of the original area A by the equation 
A* = A cos yg > or taking into account shear deformation 
ds* both in the plane and perpendicular to the plane of the 
ring leads to 



A = A cos y 2 cos 


y 3 


(16) 


Figure 3.- Cross-sectional area A 
before and after shear deformation 
y 2 of the ring. 


K the shear strains 


y 2 


and 


^3 


are small, 


A* 




(17) 


With ^2 an d ^3 known, equations ( 15 ) and ( 17 ) can be substituted into equation (8) to 
obtain the change in volume. 

An expression for y^ and 73 can be found in the following manner: The local 
shear-stress resultant Nq# can be expressed in terms of ring shear and moment result- 
ants as 


xT Mj V£ cos e V3 sin e 

019 ~ 2n-r2 + Hr irr 


(18) 


where is the twisting moment and V2 and V3 are the shear forces in the plane 

of the ring cross section. (See fig. 1 .) The energy due to shear can be written as 
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.. n2i r n2ir 

U s = g J 0 J 0 N a0 e o0 R ( 1 + r sin 0 ) r da d0 


(19) 


Substituting equations ( 6 b) and (18) into equation (19) and integrating over 6 yields the 
following result: 


! f~ Wl » + coo / 

Us = 2J 0 L M l-V^ +V2 (-“3 


u 2 * - U 1 
R 


+ V 3 0)2 + 


CJ 


R da 


( 20 ) 


The coefficients of V 2 and V 3 are the desired quantities and >"3> respectively. 
Hence 

u 2 * ~ U 1 


y 2 = -0)3 + 

u 3 - 

r3 = C 02 + — 


R 


(21) 


J 


Substituting equations (15) and (17) into ( 8 ), using equations (21), and retaining up 
to quadratic terms in the displacements gives the change in volume of the ring due to 
deformation as 


AV = 



+ 


co 3 U2’ - CO3U1 - U3'co 2 


R 



( 22 ) 


With II 2 = -p AV, the variation in the potential energy resulting from internal pressure 
can be written as 


_ r>2 7r/ui'+uo cooU<,' 

™ 2 ■ -P A 6 J„ 


- 0) 3 U 1 - Ug'COg 

R 


O 2 ' 
" 21 + ( w 3 


|R da 


(23) 


Additional potential energy of prebuckling membrane stress .- In order to take into 
account buckling behavior, the change in potential energy of the prebuckling membrane 
stresses due to a uniform line load q, as they act through the nonlinear buckling displace- 
ments, must be accounted for. The prebuckling membrane stress is taken to be 


°o = 


-qR(l + c) + pA 


(24) 


This membrane stress acts through the strain given by equation ( 6 a). Hence the variation 
of this energy while it acts through the buckling displacements is given by 

r»2ir p 27 T 

6 II 3 = 6 \ \ a 0 eQ,rt(R + r sin 0 )d 0 da (25) 

J 0 J 0 
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Substituting equation (6a) into equation (25), integrating, and retaining terms up to order 
c2 results in 


-qR(l + c) + pA r 2 * 

6I1 3 = 2 6 3 0 


u l' +u 2, ( u 3? 


+ u ^ 2 


R 


R^ 


U 2’( U 2’ - 2u l)" 
R2 


R do; (26) 


Work done by the external loads .- The toroidal ring is loaded by a uniformly distrib- 
uted radial line load q acting on the outer surface of the cross section. (See fig. 4.) 

The load q is initially directed toward the center, but 
during buckling it is assumed to behave in one of the fol- 
lowing ways: (a) The load remains in the same direction 
(constant-direction loading), (b) the load remains directed 
toward the center of initial curvature (radial loading), and 
(c) the load remains normal to the deformed surface at 
every point (hydrostatic loading). 

To evaluate the work done by three loading conditions, 
it is necessary to express the local displacements of a 
point on the surface of the shell directly under the line load 
in terms of the rigid-body motions of the cross section of 
the ring. Since e = 77 / 2 directly under the line load (see 
fig. 4), the local displacements at the point of the load can 
be obtained from equations (2) as 

u 0=tt/ 2 = U 1 - rw 3 
v 0=tt/ 2 = -™1 - u 3 
w 0=7f/2 = u 2 - |(wi 2 + ^3 2 ) J 

(a) For the constant-direction loading 
(see fig. 5) the variation in the external work 
is 

6W a = -J * * 3® ’ * 5 w e=W 2 (28) 

s 

where ds * is an infinitesimal length of the 
outer surface of the ring directly under the 
external load q. This infinitesimal length 
ds* is obtained from equations (9), (10), and 
(11) by using the displacements u, v, and w 




Figure 4.- Direction of external line 
load q. 

"N 

> (27) 
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from equations (27) in place of u^, -u 3 , and U 2 , respectively. The result is (keeping 
up to first-order terms in displacements) 


^ / Ui' - Rco) 

ds** (1 + -i 


R(1 


and the variation in the external work becomes 


m 3* + u 2\ 

+ c) / 


R(1 + c)do! 


(29) 


(b) For radial loading (fig. 6) the variation in the external work is 


(30) 


6Wb = - V q ds * 6AR 
J s* 


(31) 


where ds * is given by equation (29) and AR is the radial distance of a point on the 

outer surface of the deformed ring from the corre- 
sponding point on the undeformed outer surface. 

Retaining up to quadratic terms in the displace- 
ments gives 



AR 


= w + — — 
2(R • 


+ r) 2(R + r) 


(32) 


0=tt/2 


Figure 6.- Two-dimensional representation 
for geometry of radial loading. 


Substituting this expression for AR into equation (31) 
and making use of equations (27) and (29) leads to the 
following equation for the external virtual work: 


6W b = -qR(l + c) J 


2tt 

0 


1 + 


Uj’ - Rco> 3 + u 3 


6 w • t h 2 + “t 2 ) 


1 - c + c 2 . . . 

u l^ + u 3^ 

2 1 

R 2 


R(1 + c) 

" ^»1« 3 " u 3 w l) + c 2 (w 3 2 + W! 2 ) 


>R do: (33) 


(c) For the hydrostatic loading, the variation in the external work is 

6W C = -q 6AS2 (34) 

where AO is the net area enclosed by the deformed and undeformed ring surface at its 
outer extremity. (See fig. 7. Here, the area enclosed by the undeformed ring surface at 
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its outer extremity is 7 t(R + r) 2 .) For simplicity the out-of-plane contribution to this 
work term will not be considered, so that it is only necessary to consider the two- 
dimensional deformations shown in figure 7. Hence, the resulting buckling loads for the 

hydrostatic loading will be applicable only for in-plane 
behavior. From figure 7 is seen that AS2 can be 
represented by the equation 



AS1 = ± 


where 


and 


R = (R + r) 


1 + 


Figure 7.- Two-dimensional representations 
for geometry of hydrostatic loading. Dis- 
placements are evaluated at 6 = 7t/2. 


/3 = 


p2irf^ 2 / dfi\| . 

io L R l 1 * sir 

(35) 

1 W u 2 1 

(36) 

R + r '2(R + r)2 ; 

J 0 =ir /2 



(37) 


R + r + w 


0=1 r/2 


Substitution of equations (36) and (37) into equation (35) results in 

■•277 r 


Af2 = l(R + r) 2 f *[-1- 
2 Jq [R + 


— (u' + 2w) + - 

r (R + r) 2 


(u 2 + w2 + u'w - uw') 


da (38) 


0=ir/ 2 


After substitution of equation (38) into equation (34) and making use of equations (27), the 
variation in the work of the hydrostatic line load q is 



CWg' - CCJj2 - 



1 - c + c^ 


R^ 


-(ui 2 + u 2 2 


R 2 < 


:2 c 3 2 - 


2RCUlt0g + UJ/U2 “ U^U 2 ' - Rcu 2 0jg' + Rcu 2 *0)3j 


R da 


(39) 


Buckling equations .- The buckling equations are now obtained by substituting equa- 
tions (7), (23), (26), (30), (33), and (39) into equation (1), taking the indicated variation, 
and integrating by parts where necessary. Equating the coefficients of the arbitrary 
variations Suj, 6u 2 , 6003, Scoj, 8u> 2 , and 6U3 separately to zero yields 


-(T + l)(ui" + U2 ') + (P + S)(Ro >3 + ui - u 2 ') - Ra>3 M - q(l + c)(ui - U2 f ) 

+ q(uj - CRW3) + q^uj - u 2 ’ - cRwg) = 0 (40a) 
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(T + l)(ui f + U2 ) + (P + S)(Ro; 3' + uj_ T - 112") + Ra>3' - q(l + c^uj/ - U2”) 


+ q(ui* + U2 - cRcog’) + q(ui T + U2 - cRo^’) + q(ui’ - cRa>3 T + U2) = 0 (40b) 


-RW3” - (uj" + U2 t ) - (P + S)(-Ra>3 - ui + U2 f ) - qcRa>3 - qc(u^ + RW3) 


+ qc(-Rcog - uj + U2 r ) = 0 


(40c) 

- R ("2’ - “0 - R 4>1" + 1 - 2r)‘ 

- qcRwj - qc(-U3 + Ra>i) - qcRa>i = 0 

(40d) 

-r(u> 2” - + (p + s + fr^ 

a 3 ' + R«, 2 ) + Rr^-^j=0 

(40e) 

-(p + s + §r)(u 3 ” + R(0 2 ') + + 

u>2') + q(l + c)u3" + q(u 3 + Rccoj) = 0 

(40f) 

where 



a _ qR 3 q _ AGR 2 r _ PAR 2 

4 El 2EI " El 

T _ EAR2 JG 

El El 

(41) 


The single -underlined terms in equations (40) apply for a constant-direction load (case a), 
the double-underlined terms for a radial load (case b), and the triple -underlined terms for 
a hydrostatic load (case c). Note that the in-plane buckling behavior is not coupled with 
out-of-plane buckling. The first three equations, (40a), (40b), and (40c), describe in-plane 
buckling in terms of the in-plane variables Uj, U2, and ajg. The second three equations, 
(40d), (40e), and (40f), describe out-of-plane buckling in terms of the variables Ug, coj, 
and u>2. 


Solutions of the Buckling Equations 
A solution to equations (40) may be taken in the form 


uj = Uj sin not 

(42a) 

u 2 1 ^2 cos no; 

(42b) 

<i ) 3 = ZUg sin n a 

(42c) 

Wj = ZJj sin n a 

(42d) 
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o >2 = u >2 cos na 
u 3 = ^3 s ^ n nQ! 


(42e) 

(42f) 


where n represents the number of buckle waves and uj and Z are constants. 

In-plan e buckling .- The solution for the in-plane buckling for each of the three 
loading conditions cited previously is obtained by substituting equations (42a), (42b), and 
(42c) into equations (40a), (40b), and (40c). This results in a matrix equation of the fol- 
lowing form : 


where 



(43) 


[K] = 


(1 + T)n2 + P + S 
n(T + S + 1) + nP 
P + S + n2 


n(T + S + 1) + nP R(P + S) + Rn2 

T + 1 + n2(P + S) Rn(P + S) + Rn 

n(P + S) + n R(P + S) + Rn2 


[Z] = 


U 1 

*2 

"31 


and where [q] takes one of the following forms : 


(44a) 


(44b) 


-(1 + c) 



0 


-n(l + c) 0 

-n2(l + c) + 1 -nRc 

0 -Rc 


(44c) 



-c 

-n(l + c) 

-Rc 

®b - q 

-nc 

-n2(l + c) + 1 

-nRc 


-c 

0 

-Rc 


-c 

-nc 

-Rc 

C Q ] c = q 

-nc 

-n2(l + c) + 1 

-nRc 


-c 

-nc 

-Rc 


(44d) 


(44e) 
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The critical load for each of the three loading conditions is obtained by setting the 
determinant 


[K>[Q] 


= 0 


( 45 ) 


A first-order approximation to these critical loads can be obtained by expanding the 
determinant and keeping only the terms which are linear in q and c. (The results of 
this procedure were found to give answers within 5 percent of those obtained from the 
exact solution for c = 0.1). The lowest buckling mode is for n = 2 so that the result of 
this expansion is: 


For constant-direction loading: 


q 


cr 


4 


1 


i 2 c i 4(1 + c) i 1 + 2c 
3 S + P T 


(46) 


For radial loading: 


^cr = 


4.5 


+ — 
2 


c + 


4 + 4.5c 
S + P 


For hydrostatic loading: 


^cr 


3 


1 + 


4 + 3c 

5 + P 


(47) 


(48) 


Out-of-plane buckling.- The out-of-plane buckling solutions for the two loading con- 
ditions are obtained by substituting equations (4 2d), (42e), and (42f) into equations (40d), 
(40e), and (40f). This results in equations of the form 


where 
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(50b) 


and where 


u 3 



takes on one of the following forms : 


-n2(l + c) 0 0 



0 0 


0 


0 


-Rc 


(50c) 


Setting 


-n2(l + c) + 1 0 



0 

c 


0 

0 



leads to the buckling loads. 


Rc 

0 


-Rc 


(50d) 


(51) 


As before, a first-order approximation to these critical loads can be obtained by 
expanding the determinant and keeping only the terms which are linear in q and c. 
Also, the lowest buckling mode is for n = 2 so that the result of this expansion is 

For constant-direction loading: 



9 + 


1.125r 
P + S 


4 + 5c + 3^ * 5c + 12.375 + 13.875c , + _T — ( 0 .5 + 0.625c) 

r p + s p + s 


(52) 


For radial loading: 



12 + 


1.5r 

P + S 


4 + 4c + 1 ± 4c + 12.375 + 13.5c + _JL_(o.5 + 0.5c) 
r p+s p+s 


(53) 
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Note that in equations (46), (47), (48), (52), and (53), the internal- pressure parameter P 
and shear-stiffness parameter S always appear in combination (P + S) so that the shear 
stiffness is augmented by the internal pressure contribution. Thus P + S can be con- 
sidered as an effective transverse shear stiffness. 

DISCUSSION OF RESULTS 


The results of this analysis can be compared with results obtained by other investi- 
gators if the shear and extensional stiffnesses are assumed to be infinitely large and the 
radius ratio c is neglected. For these cases equations (46), (47), (48), (52), and (53) 
become, respectively: 


^cr 

^cr 

= 4.5 

Q 

= 3 

~cr 


_ 9 

^cr 

•44 


12 

^cr 

4 + 1 


(In-plane buckling, constant-direction loading) (54a) 

(In- plane buckling, radial loading) (54b) 

(In-plane buckling, hydrostatic loading) (54c) 

(Out-of-plane buckling, constant-direction loading) (54d) 
(Out-of-plane buckling, radial loading) (54e) 


Equation (54a) agrees with results given by Ratzersdorfer (ref. 1) and equations (54d) 
and (54e) agree with Timoshenko (ref. 3). Equations (54b) and (54c) agree with results 
given by Boresi (ref. 2). Boresi also included the effect of ring thickness, which will be 
discussed later. 

The effects of finite shear stiffness with appropriate contribution of internal pres- 
sure and radius ratio c on the buckling load can be seen by referring to figures 8 and 9. 
The results shown in these figures were obtained from the approximate equations (46), 
(47), (48), (52), and (53). However, they are quite accurate for the ranges of ring param- 
eters (c, 1/T, r, 1/(S + P)) considered here. 

In figure 8 the in-plane buckling load parameter q cr (obtained from eqs. (46), (47), 
and (48)) is plotted against the shear -stiffness parameter 1/(S + P) for c = 0 and 0,1 
with the extensional stiffness parameter 1/T equal to zero. It can be seen that the crit- 
ical load is influenced significantly by changes in the direction of loading during deforma- 
tion. It is also seen that the critical load is influenced by c and that this influence is 
likewise dependent on the direction of loading. For example, as c increases q is 
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Figure 8.- Variation of in-plane critical load parameter with effective shear- 
stiffness parameter. 

decreased. The influence of ring depth on the buckling load parameter shows the same 
trends that were obtained by Boresi (ref. 2) for the radial and hydrostatic loading condi- 
tions of a ring of rectangular cross section. However, there appears to be no literature 
available for the constant -direction loading with the effects of ring depth included. Cal- 
culations were also carried out for finite but realistic values of 1/T in the range 
0 S s 0.01 for the constant-direction loading case (eq. (46)). Throughout this range, the 
values of q cr changed so slightly that the differences could not be plotted in figure 8. 
Thus, ring extension has little influence on in- plane buckling in the n = 2 mode. This 
would not be the case, however, if external constraints forced buckling in a higher mode. 

Note that in all cases, the critical load parameter markedly decreases with an 
increase in the shear stiffness parameter l/S + P. Furthermore, since the pressure 
parameter P always appears with the shear parameter, the internal pressure will have 
no effect if the shear stiffness of the membrane wall is large (S = °°). 

In figure 9, the out-of-plane buckling load parameter q cr (obtained from eqs. (52) 
and (53)) is plotted against shear stiffness parameter 1/S + P for c = 0 and 0.1 with 
twisting stiffness r = 1.0. It can be seen that, as in the in-plane buckling case, the 
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Figure 9.- Variation of out-of-plane critical load parameter with effective 
shear-stiffness parameter for r = 1. 

direction of loading during deformation, as well as the value of c, significantly affects 
the magnitude of the buckling load. In all cases the buckling load parameter decreases 
with an increase in 1/S + P. Also, the out-of-plane buckling loads are lower than the 
in- plane buckling loads for the same external loading conditions and geometrical proper- 
ties. Note in equations (54d) and (54e) the potentially large effect of the parameter r on 
out-of-plane buckling. Evidently, if r = 1, drastic reductions in the buckling load occur. 
Such low values may be possible for deployable toroidal rings constructed from single 
layers of organic or metallic fabric (see, for example, ref. 6). 

It is interesting to note that the in- plane and out- of- plane buckling modes are uncou- 
pled. This is true because the product of inertia for the toroidal ring is zero. A conse- 
qence of this uncoupling is that the twisting stiffness T does not affect in-plane buckling 
and the extensional stiffness T does not influence out-of-plane buckling. 

CONCLUDING REMARKS 

A theoretical investigation of the buckling of a pressurized toroidal ring under a 
uniformly distributed line load has been carried out, with finite shear stiffness, 
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extensional stiffness, and the appropriate contributions from the internal pressure taken 
into account. The equations governing both in-plane and out-of-plane buckling behavior 
of the ring for three directions of loading have been derived and solved, and simple buck- 
ling formulas are given. 

The major results of the investigation are: 

1. The magnitude of the buckling load is very sensitive to the direction of loading 
during the buckling deformations, and also to the ratio of cross-sectional radius to 
radius of the ring. 

2. Out-of-plane buckling loads are lower than the in-plane buckling loads for the 
same external loading conditions. 

3. Shear stiffness of the ring wall material and internal pressure both contribute to 
the effective transverse shear stiffness of the ring, which strongly influences buckling 
strength. 

4. The effect of extensional stiffness on the in- plane buckling load is negligible. 
However, twisting may strongly influence the buckling loads. Out-of-plane buckling 
behavior does not couple with in-plane buckling behavior. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., March 3, 1967, 

124-08-06-21-23. 
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